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Abstract 

This paper reports the results of Suzaku observation of the speetral variation of the 
blaek hole binary LMC X-1 in the soft state. The Observation was earried out in 2009 
from July 21 to 24. The obtained net eount rate was ~ 30 eounts s“^ in the 0.5-50 keV 
band with ~ 10% peak-to-peak flux variation. The time-averaged X-ray speetrum eannot 
be deseribed by a multi-eolor disk and single Compton eomponent with its refleetion, but 
requires an additional Comptonized emissions. This double Compton eomponent model 
allows a slightly larger inner radius of the multi-eolor disk, implying lower spin parameter. 
Signifieant speetral evolution was observed above 8 keV along with a flux deerease on a 
timeseale of ~ 10^-10^ s. By speetral fitting, we show that this behavior is well explained 
by ehanges in the hard Comptonized emission eomponent in eontrast to the maintained disk 
and soft Comptonized emission. 

Key words: aeeretion, aeeretion disks — blaek hole physies — X-rays: individual 
(LMC X-1) 


1. Introduction 

X-ray spectra from stellar mass black hole (BH) binaries are known to transit between the 
high/soft state (HSS) and low/hard state (LHS), in accordance with the mass accretion rate. The 
HSS spectrum is mainly composed of a bright thermal component that dominates the soft X-ray 
band accompanied by a less bight power law (PL) component that exhibits a “hard tail” in the ^10 
keV band. The thermal component is well reproduced by a multi-color disk (MCD; Mitsuda et al. 
1984; Makishima et al. 1986) emission from an optically thick and geometrically thin accretion disk 
described by the standard disk model (Shakura & Sunyaev 1973). From the fact that the inner radius 
of the accretion disk derived from the MCD model is constant (e.g., Ebisawa et al. 1993; Steiner et 
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al. 2010), the existence of an innermost stable circular orbit (ISCO), which is determined by the BH 
mass and spin, has been established (see also Makishima et al. 2000). 

The hard-tail component seen in the HSS is thought to be produced by Compton upscatter- 
ing of the MCD photons by energetic electrons. This implies a situation similar to the one that 
produces the dominant hard X-ray continuum in the LHS, which is understood as resulting from ther¬ 
mal Comptonization of the disk photons by hot thermal electron clouds (e.g., Gierlinski et al. 1997; 
Frontera et al. 2001; Zdziarski & Gierlinski 2004; Makishima et al. 2008; Takahashi et al. 2008). 
However, these two phenomena show considerable differences. The hard tail in the HSS is expressed 
by a single PL with a photon index of T ~ 2.0-2.5, and often extends to ~ 1 MeV, in contrast to 
the LHS continuum which is flatter (T < 2) with clear a cut-off at ~ 100 keV. As a result, the HSS 
hard tail can be explained either by non-thermal Comptonization, or by a hybrid of non-thermal and 
thermal Comptonization processes (Zdziarski et al. 1993; Gierlinski et al. 1997; Coppi& Maccarone 
2002). Overall, the hard-tail phenomenon is much less well understood than the soft component. 

When a BHB (particularly a transient one) evolves from the LHS into the HSS, it often passes 
through another spectral state, called the very high state (VHS; e.g, Miyamoto et al. 1991). The VHS 
spectrum is characterized by luminous disk emissions with a strong PL hard tail with a steep slope of 
r > 2.4. The latter is resolved to two hard-tail, and explained by non-thermal plus thermal combined 
Comptonization (e.g., Kubota & Done 2004 ; Kobayashi et al. 2003), of which the thermal fraction 
is thought to decrease as the system approaches the HSS from the VHS. This state evolution is also 
thought to involve changes in the optically-thick disk, which extends down to ISCO in the HSS (as 
described above) but while is likely to be truncated in the LHS (e.g., Makishima et al. 2008). The 
VHS is thus understood as a transient stage in which the innermost disk radius gradually decreases 
and approaches the ISCO. (Kubota & Done 2004; Done & Kubota 2006; Tamura et al. 2012). In this 
way, the optically-thick disk and the Comptonizing clouds are suspected to be coupled to each other 
through the state transition. In the present paper, we examine possible relationships between the disk 
and the Compton-clouds in BHBs by examining a VHS that is rather close to the HSS. This objective 
requires broadband coverage from ~ 1 to ~ 100 keV, for which Suzaku is ideal. 

LMC X-1 is a persistently X-ray luminous BH binary that is accompanied by an O type star 
with an estimated mass of 32 Mq, where Mq is the solar mass. The binary parameters are known 
with good accuracy: BH mass Mbh = 10.9 Mq, inclination angle i = 36°4, and distance 79 = 48 kpc 
(Grosz et al. 2009). For decades, it has been found in the HSS or VHS, with relatively high luminosity 
of ~ 10% of the Eddington value. Previous observations have reported a variable inner disk radius and 
a steep hard-tail slope of F ~ 3, which are both suggestive of the VHS. This, together with its strong 
X-ray variations (e.g., Ruhlen et al. 2011), makes the object a good source for our study. In previous 
studies of LMC X-1 in energy bands below 20 keV, the hard tail was in most cases reproduced by a 
simple PL feature with reflection (e.g. Gou et al. 2009; Ruhlen et al. 2011). 

Steiner et al. (2012) showed clear evidence for an iron line in the time-averaged Suzaku spec¬ 
trum, and a positive correlation between the iron line flux and Compton scattering fraction obtained 
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from RXTE observations. Although those results suggest that the Comptonized hard-tail photons 
are illuminating the disk to produee the iron line photons, that study considered neither combined 
Comptonization nor studied spectral variations in the Suzaku data. In the present paper, we hence 
reanalyze the same Suzaku observation data, attempting to see whether the inner disk radius remained 
constant, and whether the hard tail exhibited the signature of combined Comptonization. 

2. Observation and Data Reduction 

Suzaku, which is the fifth Japanese X-ray satellite (Mitsuda et al. 2007), carries onboard the 
X-ray imaging spectrometer (XIS; Koyama et al. 2007) and the hard X-ray detector (HXD; Takahashi 
et al. 2007; Kokubun et al. 2007). The XIS consists of four charge-coupled device cameras placed 
at the focus of the X-Ray Telescope (XRT; Serlemitsos et al. 2007), covering the 0.2-12 keV energy 
range. However, one of four units is not operational (SUZAKU-MEMO 2007-08'). Of the three 
available XIS units, two (XIS-0 and XIS-3) are front-illuminated (El) while the other (XIS-1) is back- 
illuminated (BI). The HXD consists of PIN silicon diodes (HXD-PIN; 10-70 keV) and Gd 2 Si 05 Ce 
scintillators (GSO; 50-600 keV). In 2009 from July 21 UT 18:38 to July 24 21:29, Suzaku observed 
the source at “XIS nominal” pointing position. The XIS was operated in standard clocking mode and 
“1/4 window” option in order to attain a time resolution of 2 s. 

We used data product from Suzaku pipeline processing version 2.4.12.26 with calibration ver¬ 
sion hxd20090511, xis20090402 and xrt20080709, and software version HEADAS 6.6.2. XIS and 
HXD events were screened by standard criteria. We discarded events collected with Earth elevation 
below 5°, or with the sun irradiated-Earth elevation below 20° (for XIS), or when the spacecraft was 
in an orbit phase within 436 s after (for XIS) or 180 s before and 500 s after (for HXD) the South 
Atlantic Anomaly ingress/egress, or at low-cutoff rigidity regions below 6 GV. We accepted only the 
XIS events with standard grades (0, 2, 3, 4 and 6) in the analysis. 

XIS spectra and light curves were extracted from a rectangle region of 8(6x4f5, attained in 
the 1/4 window area and centered on the image peak, while the background was also obtained from 
a source-free region in the same 1/4 window area. Eurthermore, we excluded a circular region of 50" 
around the image peak. According to the software AEPILEUPCHECKUP (Yamada et al. 2012), the 
remaining signals are evaluated to have a pileup fraction of <1%. The net averaged XISs count rates 
were ~15 counts s“l in the exposure time of ~ 110 ks for each sensor. The XIS redistribution matrix 
files were calculated using the xisrmfgen tool and the ancillary response files were simulated by 
the xissimarf gen tool (Ishisaki et al. 2007). 

The XISl unfolded data deviate systematically, at energies of ^ 3 keV, from those of XISO 
and XIS3. This is presumably due to the residual calibration uncertainties of XISl (Ishida et al. 2011, 
subsection 7.3.3 of Suzaku technical description^ ). Therefore we did not use the XISl data in the 

' Available at http://www.astro.isas.ac.jp/suzaku/doc/suzakumemo/suzakumemo-2007-08.pdf 
^ Available at http://www.astro.isas.ac.jp/suzaku/doc/suzaku_td/ 
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following analysis. 

The HXD screened events were used to obtain spectra and light curves from the PIN and 
GSO sensors. The net exposure of each sensor was 124 ks after dead-time correction. The dead¬ 
time fractions were 6.8 % averaged for both PIN and GSO sensors. The cosmic X-ray background 
was modeled assuming an exponentially cutoff PL model (Boldt et al. 1987). Non-X-ray back¬ 
ground (NXB) models were provided by the HXD team (Fukazawa et al. 2009). We used the 
model with METHOD=“LCFITDT(bgd_d)” , together with the version of METHODV=“2.0ver0804” 
and “2.4ver0912-64” for the data from the PIN and GSO sensors, respectively. Eigure 1 shows the 
background-subtracted spectra. Events are significantly detected by the PIN sensor up to 50 keV, 
which is above the systematic uncertainty in the NXB model (3% for PIN). 

In addition to EMC X-I, there were three hard X-ray sources in the PIN field of view (EOV; 
30'x30'); PSR B0540-69.3, SN 1987A and RX J0541.4-6936. Although the expected 20-100 keV 
flux of PSR B0540—69.3, ~ 2.9 x 10“^^erg cm“^ (Campana et al. 2008), is comparable to the 
observed PIN flux of ~ 5.5 x 10“^^erg cm“^ (modeled with a single PE), it is separated from 
EMC X-1 by 25'. The PIN EOV during the present observation is evaluated to have a transmission 
efficiency of 21%. Therefore, in the following analysis, the spectral contribution from this active 
pulsar was modeled by a PE with r=2.12 (Campana et al. 2008), but with the flux reduced to 21% 
of that reported. In the light curve analysis, we subtracted the constant count rate of 1.21x10“^ 
counts s“^ in 13-50 keV, which estimated by PIMMS^ tool, from the PIN light curve. The hard X-ray 
spectra of SN 1987A and RX J0541.4—6936 were cutoff-PE with r=1.7, cut-off energy of i7c=100 
keV and flux of ~ 10“^^ erg cm“^, and r=1.0, i7c=15 keV and the flux of ~ 10“^^erg cm“^, 

respectively, as taken from the INTEGRAL source catalog ^ (Bird et al. 2010). Their locations are 
off-axis to the PIN EOV, and the PIN efficiency was <50% of the peak value. The contributions of 
these two sources are thus an order of magnitude below the PIN flux and hence negligible. Because 
the GSO EOV contained too many contaminating sources to individually evaluate, we do not use the 
GSO data in this paper. 

3. Analysis of Time-Averaged Spectra 

Using the X-ray spectral fitting package XSPEC version 12.8.2 and employing the solar abun¬ 
dance tabulated in Anders & Grevesse (1989), we analyzed 0.8-10 keV XIS spectra and 13-50 keV 
HXD-PIN spectra, averaged over the entire exposure. The energy bands of 1.6-2.0 keV and 2.2-2.4 
keV were excluded to avoid relatively large systematic instrumental uncertainties near the Silicon K 
edge (1.74 keV and 1.84 keV for Ka and K/3, respectively) and the Gold M edge (2.29 keV), re¬ 
spectively. The remaining instrumental feature at ~ 3.2 keV by the Gold M edge was modeled by a 
Gaussian line with a fixed width of a = 0.1 keV (ref. Kubota et al. 2007). We used a combination of 

^ Avaliable at http://heasarc.gsfc.nasa.gov/docs/software/tools/pimms.html 
^ Available at http://isdc.unige.ch/integral/science/catalogue 




Fig. 1. Background-subtracted spectra of LMC X-1 obtained with Suzaku XISO (black), XIS3 (red) and HXD-PIN 
(green). The blue line indicates the systematic uncertainty level of PIN. 

XISO and XIS3 after co-adding them together, and we added a 1% systematic uncertainty to all XIS 
energy channels. We extended the energy range of spectral fitting to 0.1-1000 keV, as used in some 
convolution models. We fitted the XIS and PIN spectra simultaneously with a cross normalization 
factor of 1.16 between XIS and PIN (SUZAKU-MEMO-2008-06^). 

We assumed the situation where part of the MCD photons are Compton-upscattered to the 
hard tail, and with reflection of Comptonized emissions. We utilize a convolution model simpl 
which takes some fraction of the disk seed photons and upscatters these to higher energy (Steiner 
et al. 2009). In order to account for the reflection of the Comptonization continuum from ionized 
material, we employed the rfxconv model (Kolehmainen et al. 2011) which is a convolution 
model and calculates reflection continuum and fluorescence based on the tables of Ross & Fabian 
(2005). We fixed the ionized parameter ^ =1000 erg cm s“^, inclination angle of 6*=36.38°(Orosz 
et al. 2009) and assumed solar abundances. Using kdblur (Laor 1991) the reflected spectrum 
is smeared by general relativistic effects, with fixed inner and outer disk radii of = 2 Rg 

(Gou et al. 2009) and = 400 Rg, inclination angle i = 36.4°, and emissivity index f3 = —3 

(Rg = GM/c^, where M is the mass of the central BH). The output of simpl*diskbb is the MCD 
component reduced by Comptonization and upscattered emissions. To calculate the reflection of the 
Compton component, we divided the simpl into simplsouj-ce and sirnplcompton for reduced MCD 
and Comptonization emission, respectively. Consequently, we employed the continuum model as ex¬ 
pressed by (simplsource + kdblur*rf xconv*simplcompton)*diskbb modified by photoelectric 
absorption (phabs in XSPEC) as Model 1. 

Figure 2 (a) shows the folded spectra and residuals from the best-fit model. The best-fit pa¬ 
rameters are shown in table 2 as Model 1. This fit is unacceptable, with x^/d.o.f = 290.13/118, and 
we see large residuals below 3 keV and above 20 keV. To check that the assumption of ^ = 1000 erg 

^ Available at http://www.astro.isas.jaxa.jp/suzaku/doc/suzakumemo/suzakumemo-2008-06.pdf 
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cm s“^ does not affect the result, we set ^ to be free. This did not resolve the problem, with results 
of x^/d-O-f = 216.84/117, where ^ x 10“^ erg cm s“^ and reflection strength of f2/27r = 

6.07lo;io, which is too large to be physically compatible. 

To introduce a physical explanation for the excess in the 5-10 keV band, we replaced the 
diskbb model by bhspec model (Davis et al. 2005) (Model 2) to account for relativity. Figure 2 
(b) shows the folded spectra and residuals from the best-fit model. The best-fit parameters are shown 
in table 2 as Model 2. The fit is still unacceptable with x^/d.o.f = 302.24/118. As shown by the black 
and red lines in figure 2 (d), the application of relativistic effects makes the disk component slightly 
broader but the difference is compensated for by increased photoelectric absorption. We tested the 
fit again by allowing ^ to be free, and derived x^/d-o.f = 245.15/117, where ^ erg cm 

s“^and f2/27r = 8.89l^o5^^- 

To provide an alternative picture, we introduce a combined hard and soft Comptonization 
model. For the soft Comptonization process, we employed a dkbbfth model (Done & Kubota 
2006), which calculates Comptonized MCD spectrum from the accretion disk and the thermal plasma 
covering the inner region at Rm <R< i?tran, where R^ and i?tran are inner disk radius and outer radius 
of the thermal Comptonizing cloud. The model parameters are temperature of seed MCD, electron 
temperature kT^, photon index Fth, 7(tran> and normalization. We again include the simpl model for 
hard Comptonization with fixed F = 2.1. We assume that hard Comptonization occurs independently 
of soft Comptonization (i.e., hard electrons do not scatter soft Comptonized photons) and assume 
the same reflection process for both components. This thus gives a continuum model expressed 
by phabs * [simplsource * dkbbfthdisk + kdblur * rfxconv * (simplnc * dkbbfthdisk 
-I- dkbbfthsc)] (Model 3), where the subscripts of HC and SC denote the hard and soft Compton 
component, respectively. 

The fitting result with the third model is shown in figure 2 (c) and table 2 as Model 3. The 
reduced is acceptable with 111.57/116. Equation 1 of Tamura et al. (2012) gives the relationship 
between kT^, Tth and optical depth assuming a source of slab geometry with seed photons at the 
bottom of the slab: 
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where 0 = kTe/mc^ and the optical depth is approximately half that of spherical geometry (see also 
Zdziarski et al. 1996). Using this equation, the optical depth is calculated to be r = 0.67^oSince 
the normalization of dkbbfth is not equivalent to that of diskbb, we calculated the inner disk 
radius pn via equation A.l in Kubota & Makishima 2004 using the unabsorbed photon flux. We 
derived r-m = 55.01J g which is consistent with the estimation from Model 1. This implies that the soft 
Compton cloud localized within /(trans = 8.9i7in scatters 12% of disk photons. 

As show in figure 2 (d), the inclusion of the soft Comptonized emission shifts the disk com¬ 
ponent to lower temperature, which makes derived larger substantially. The fraction of hard 
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Comptonized emission decreases and its photon index flatten, due to the soft Compton component. 


Table 1. Best-fit parameters for time averaged spectra. 




Model 1 

Model 2 

Model 3 

phabs 

Vh (xlO^i cm2) 

r ^y+0.07 

O.^/0.06 

^ qo+0.07 

O.UO o.o6 

6.06///i 

diskbb 

fcTin (keV) 

0.784 ±0.006 

— 

— 


Norm * 

104.69/3-94 

— 

— 

bhspec 


— 

0.79 ±0.01 

— 


L/LEdd 

— 

0.137±0.001 

— 

dkbbfth 

Tin (keV) 

— 

— 

noil +0.012 

U.»li_o 010 


Norm (X10-2) 

— 

— 

1 QQ+O.OO 
-^■'^^-0.12 


-^tran (-^in) 

— 

— 

8.86//°® 


fcTe (keV) 

— 

— 

15.5//I 


Tth 

— 

— 

q Q1+011 
'^•^-'--0.10 

rfxconv ^ 

tl/27r§ 

0.64/°:“ 

0.64/°:“ 

Q QS+0-34 

simpl 

/scat '' 

0 1 (37+0022 
U.iO/_Q 015 

0.164/°:°33 

0 nQo+0.002 

U.Ud)d)_o.o07 


r 

^■°®-0.05 

2 QO+0 04 
^•yu o,o5 

2.1 (fixed) 

xVd.o.f 


290.13/118 

302.24/118 

111.57/116 

Derived value 

Tin (km) * 

54.8//^ 

47.6 ±0.8 

55.0/1:° 


Tth 

— 

— 

0.67/°:°° 


/th 

— 

— 

0.12 


Notes. The errors are 90% confidence level for single parameter. 

* Normalization defined as r?^cos2/(D/10 kpc)^ 

^ The solar abundances assumed, with fixed ionization parameter of ^ = 1000 
^ Reflection strength of the iref lect model. Q is solid angle of the reflector. 

II Scattering fraction of the simpl model 
^ Inner disk radius derived from fit results. 

Model 1: Calculated from the normalization of diskbb. 

Model 2: Equivalent to a* of bhspec calculated with Mbh = 10.9Mq. 

Model 3: Estimated from unabsorbed photon flux via Equation (A.l) in Kubota & Makishima (2004). 


4. Light Curves and Spectral Variation 

Figure 3 shows background-subtracted light curves and hardness ratio variations for LMC X- 
1. Both the X-ray counts and hardness ratios exhibit variations on a time scale of 10^-10^ s. In 
order to estimate the energy dependence of the variation, we calculate the rms value for each divided 
energy bands using a bin size of 5760 s, which is the timescale shown in figure 3. The derived rms 
spectrum is shown in figure 4. The fractional rms variation clearly depends on the energy band, 
increasing towards higher energy range and exceeding 10% above ~6 keV. 
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Fig. 2. (a) Top: Folded spectra of LMC X-1 from the best-fit for phabs * [ ( simplsource + kdblur * rf xconv 
* simplCompton ) * diskbb ] (Model 1). Black and red plots (and lines) denote the spectra (and best-fit models) 
for XIS(0 h- 3) and HXD-PIN, respectively. Bottom: Residuals between the data and best-fit models. 

(b) The same as (a), but for phabs * [ ( simplsource + kdblur * rf xconv * sirnplcompton ) * bhspec ] 
(Model 2). 

(c) The same as (a), but for phabs * [ simplsource * dkbbfthaisk + kdblur * rfxconv * ( simplnc * 
dkbbfthdisk + dkbbfthsc ) ] (Model 3). 

(d) Unabsorbed model spectra from the best-fit Model 1 (black), Model 2, (red) and Model 3 (blue). Emission 

from the disk (simplsource*diskbb, simplsource*bhspec and simplsource*dkbbfthdisk). Compton 
(simplcompton*diskbb, simplcompton*bhspec, and simplHC*dkbbfthdisk) and reflection components 
are represented by dotted, dash-dotted, and dashed lines, respectively. For Model 3, the soft Compton component 
(dkbbfthsc) and the reflection of that are represented by light blue. 


Figure 5 shows the hardness ratio of 6-10 keV to 4-6 keV against 4-10 keV eount rate with 
a time bin size of 5760 s. The hardness ratio is elearly eorrelated with souree flux, whieh indieates 
signifieant speetral variation. In order to evaluate the speetral variation shown in figure 3 and 4, we 
divided the data into hard phase (HP) and soft phase (SP) for hardness ratios higher and lower than the 
average value of 0.25 (dashed line in figure 5). Figure 6 shows speetra normalized by PL with F = 2 
for eaeh phase and pulse-height ratio of HP/SP The differential speetrum has a peak at around 6 keV 
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Fig. 3. Background-subtracted light curves of the XIS FI in the 0.5-10 keV band and of the HXD-PIN in the 13-50 
keV band, and hardness ratio variation of PIN/XIS. The bin size is set to a Suzaku orbital period of 5760 s. The 
constant flux from PSR B0540—69.3 is subtracted from PIN light curve. 



Fig. 4. Fractional rms variation spectra for XIS FI and HXD-PIN. The bin size is set to 5760 s. In the calculation 
of data points of 13-20 keV and 20-50 keV, the constant flux from PSR B0540—69.3 is subtracted from PIN light 
curve. 

and dominates above~10 keV. As the pulse-height ratio shows, although the variation is < 10% in 
the band below 4 keV, it is ~ 30-50% in the band above ~8 keV. In the time-averaged speetra (figure 
2 ), the eontribution of Compton and refleetion eomponents dominates above ~6 keV in the energy 
speetrum with a peak around ~4-6 keV. The behavior thus implies speetral variability of Compton 
and/or reflection components. 

To quantify the variations in the continuum spectra, we performed a broadband spectral fit¬ 
ting to each phase using Model 3. The best-fit parameters are shown in table 2. The fit for both 
HP and SP are acceptable, with x^/d-o.f = 109.95/116 and 108.64/116 for HP and SP, respectively. 
The hard Compton scattering fraction /scat differs significantly by ~40% from O.OdSlg ogy/HP) to 
0.027lo oog/SP). In contrast, the MCD component and photo-electric absorption were stable between 
two phases. This suggests that the hard Compton cloud changes independently from the disk condi- 
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Fig. 5. Hardness ratio of 6-10 keV to 4—6 keV against 4-10 keV count rate for a bin size of 5760 s. The dashed 
line represents the average value of hardness ratio of 2.5 and indicates the limit between the hard and soft phase. 


tion. Although the soft Comptonization fraction fth decreases slightly from 0.12 (HP) to 0.11 (SP), 
the changes in the other parameters i?tran, kTg, Fth and rth are insignificant. 


(a) 



(b) 



Fig. 6. (a) Spectra divided with PL (L = 2.0) of hard phase (black), soft phase (red) and difference between the two 
phases (green), (b) Ratio of the spectra of the hard phase to the soft phase. 


5. Summary and Discussion 

We analyzed Suzaku data of LMC X-1 observed in 2009 from July 21 to 24 and obtained 0.8- 
50 keV spectra and light curves. The time-averaged spectra cannot be explained by simple modeling 
of a non-relativistic or relativistic disk and single Comptonized emission (Models 1 and 2 in table 2) 
and the spectra exhibit large deviations below 3 keV, around 6 keV, and above 20 keV. These residuals 
are successfully compensated for by adding a hard and soft Comptonization model (Model 3). The 
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Fig. 7. Unfolded spectra from the best-fit model (Model 3) for the (a) hard phase and (b) soft phase. 

(c) Unabsorbed model spectra of HP (blue) and SP (red) from the best fit to Model 3. Emission from the disk 
(siinplsource*dkbbfthdisk). Compton (siinplcompton*dkbbfthdisk) and reflection components are repre¬ 
sented by dotted, dash-dotted, and dashed lines, respectively. The soft Compton compo-net (dkbbfthcompton) 
and its reflection are represented by light blue and orange for HP and SP, respectively. 


hard and soft components are deseribed by a PL tail with F = 2.1 representing the exeess above 20 keV 
from the simpl model and thermal Comptonization at inner disk region of ^ 9 Rm with kT^ ~ 16 
keV and rth ~ 0.7 to deseribe the 4-10 keV residuals by the dkbbf th model (figure 2), respeetively. 
Using the best-fit parameters, the unabsorbed 0.8-50 keV flux is estimated to be 8.7 x 10“^° ergs em“^ 
s“^, whieh gives an X-ray luminosity of 2.4x 10^® ergs s“^ for an isotropie emission from 48 kpe. The 
bolometrie luminosity was also ealeulated to be 3.3 x 10^® ergs s“^. 

In eontrast to the ease where the speetrum is not fitted by the single Compton model, Steiner 
et al. (2012) reported redueed by similar fitting to Model 2 from the time-averaged Suzaku 

speetrum of LMC X-1. This differenee eould be mainly eaused by absorption model for inter-stellar 
medium (ISM). They employed the absorption model of TBvarabs (Wilms et al. 2000) with ISM 
eomposition taken from Hanke et al. (2010), while we employed phabs in the analysis. Replaeing 
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Table 2. Best-fit parameters for each phase. 




HP 

SP 

phabs 

TVh (x 10^^ cm^) 

e.OQtna 

6.01 ±0.08 

dkbbfth 

kTin (keV) 

0.813tni2 

0.807t°;°°I 


Norm (x 10“^) 

1 qr+0.10 

1.42lHo 


-^tran (-^in) 


Q 07 + 1-02 
-0.92 


/cTe (keV) 

M.Olt.? 

16.31+® 


Fth 

'^•^^-0.14 

^■^^-0.26 

rfxconv ^ 

n/27r§ 

0.85 ±0.24 

u.oy o.36 

simpl 

^scat 

o.o48;°:r7 

u.uz ( _o.o09 


r 

2.1 (fixed) 

2.1 (fixed) 

x'^/d-o.f 


109.95/116 

108.64/116 

Derived value 

Tin (km) * 




Tth 

0.70 ±0.16 

0-6i;°;?^ 


fth 

0.12 

0.11 


Notes. The errors are 90% confidence level for single parameter. 

^ The solar abundances and ionization pai'ameter of ^ = 1000 assumed. 

^ Reflection strength of the rf xconv model. 17 is solid angle of the reflector. 

II Scattering fraction of the simpl model 

^ Inner disk radius estimated from unabsorbed photon flux via Equitation (A.l) in Kubota & Makishima (2004). 

the phabs to TBvarabs, we obtain improvement of redueed xl from 2.6 to 1.5. However, Model 
3 still gives the best fit results with the redueed xl of 1-1 • Therefore these differenee of absorption 
models dose not affeet to the conelusion of this paper. 

We found signifieant variation both in the XIS and PIN band. The eount rates eorrelate with 
the hardness ratio, and different spectral changes appeared in the three energy bands of below 2 keV, 
above 10 keV, and 2-10 keV, which exhibited variations of ~1% , ~40%, and values in between, 
respectively (figure 6b). This implies the existence of three spectral components, which is consistent 
with the modeling of MCD and combined hard and soft Comptonized emission. From the spectral fit¬ 
tings, we found that the spectral change is mainly explained by variations in the hard Comptonization 
fraction by a factor of ~ 40% in contrast to keeping the MCD parameters of Tin and Rin- 

The obtained results show that the spectrum is composed by three spectral components of disk 
and two Comptonized emissions instead of disk with single Comptonized emission. This increases 
inner disk radius from 47.6±0.8 km (Model 2) to 55.01^} g km (Model 3), hence, implies lower BH 
spin. 

Figure 8 shows a comparison of the best-fit parameters of soft Comptonized emission from 
three VHS BHBs observed by Suzaku, GX 399—4 (Tamura et al. 2012), 4U1630—47 (Hori et al. 
(2014)), and LMC X-1 (this work). The parameters of kT^ and fth derived from LMC X-1 are 
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lower than those of the other studies, although rth is similar. We also found a stable rin at ~55 km 
(corresponds to ~3-4 Rg with M = 10.9 Mq), which is assumed to reach the ISCO, of LMC X-1. In 
contrast, the truncated inner disk radii of GX339—4 and 4U1630—47 were reported by Tamura et al. 
(2012) and Hori et al. (2014), respectively. Taking into account that the inner disk radius and the soft 
Comptonized emission decrease through the state evolution from VHS to HSS, we found the feature 
that LMC X-1 is in the VHS rather than HSS, and the inner disk radius reaches the ISCO though the 
relatively weak soft Comptonized emission. 



0.1 0.15 0.2 0.25 0.3 


7h 


Fig. 8. Comparison of derived dkbbfth parameters of rth, kTg and fth from three objects. Squares, circles and 
asterisks represent the parameters of GX 399—4(Tamura et al. (2012)), 4U1630—47 (Hori et al. (2014)) and LMC 
X-1 (this work), respectively. 
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